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We consider the gapless modes along the vortex line of the fully gapped, momentum independent 
paired states of three dimensional Dirac fermions. For this we require the solution of fermion zero 
modes of the corresponding two dimensional problem in the presence of a point vortex, in the plane 
perpendicular to the vortex line. Based on the spectral symmetry requirement for the existence of 
the zero mode, we identify the generalized Jackiw-Rossi Hamiltonians for different paired states. A 
four dimensional generalized Jackiw-Rossi Hamiltonian possesses spectral symmetry with respect to 
an anti-unitary operator, and gives rise to a single zero mode only for the odd vorticity, which is 
described by a Zi index. In the presence of generic perturbations such as chemical potential, Dirac 
mass and Zeeman coupling, the two dimensional problem for the odd parity topological supercon- 
ducting state maps onto two copies of generalized Jackiw-Rossi Hamiltonian, and consequently an 
odd vortex binds two Majorana fermions. In contrast there are no zero energy states for the topolog- 
ically trivial s-wave superconductor, in the presence of any chiral symmetry breaking perturbation 
in the particle-hole channel, such as Dirac mass. We show that the number of one dimensional 
dispersive modes along the vortex line is also determined by the index of the associated two dimen- 
sional problem. For an axial superfiuid state in the presence of various perturbations, we discuss 
the consequences of the Zi index on the anomaly equations. 



I. INTRODUCTION 

The existence of zero modes for the Dirac fermions in 
the presence of a topologically non-trivial configuration 
of an order parameter or a gauge field is an interesting 
problem in condensed matter and high energy physics. 
In a seminal paper, Jackiw and Rebbi demonstrated 
the emergence of the fermion zero modes on various 
defects in odd spatial dimensions^. In particular, they 
have showed the existence of zero mode for a domain 
wall in one dimension, t'Hooft-Polyakov monopole and 
dyon in three dimensions. The fermion zero modes due 
to a solitonic defect can give rise to induced quantum 
numbers and fractionalizationi^. The domain wall of a 
scalar order parameter (Dirac mass) in one dimension 
binds single zero mode, leading to the fractionalization 
of charge, which has been observed experimentally in 
polyacetylcncA Interestingly, the edge and the surface 
states of many gapped topological systems in higher 
dimensions arc also determined by one or multiple copies 
of the one dimensional Jackiw-Rcbbi model. Some 
specific examples are the edge states of the d + id A 
and the p + ip%& superconductors and the quantum 
spin Hall stated in two dimensions, and the surface 
states of the three dimensional topological insulators 
and superconductors^. For these higher dimensional 
problems, the domain wall describes a boundary between 
topologically distinct vacua. It is also interesting to note 
that the chiral surface states of gapless Weyl semi-metal 
and superconductors are also governed by similar one 
dimensional problem^. 

The point vortex and the line vortex of a U(l) sym- 
metry breaking order parameter respectively in two and 
three dimensions, arc also interesting topological defects. 



The point vortices of the p + ip superconducto r 5 ' 10 ' 11 
and the Kckule bond-density wave order in graphenei^ 
support fermion zero modes. In contrast the line 
vortex of many three dimensional paired states such as 
axial superfluicU^, axial superconductor—, superfiuid 
3 He — B^, the gapped superconducting states of three 
dimensional Dirac fermionsi^— support one dimensional 
gapless fermions along the vortex core. If we consider, 
the line vortex along the z direction, the momentum k z 
is a conserved quantity, and the dispersing states are 
localized in the x — y plane. Therefore, for k z = 0, we 
indeed solve a two dimensional problem of the fermion 
zero modes. Consequently, the number of dispersing 
modes becomes equal to the number of fermion zero 
modes of an effective two dimensional problem in the 
presence of a point vortex. 

The problem of zero energy modes bound by a point 
vortex in two dimensions has been considered by Jackiw 
and Ross i 20 ' 21 . The pertinent Hamiltonian is 
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where cfy's are the derivatives with respect to the spa- 
tial coordinates, and Ai,A2 are the components of a 
complex order parameter. The four dimensional Her- 
mitian r-matrices satisfy the anti-commuting algebra 
{I^r,-} = 2%, where i,j = 1,2,3,4. The gauge field 
corresponding to the broken C/(l) symmetry is denoted 
by a.,-, and the finite amplitude of the order parameter 
gives rise to the Meissner effect for a,j. The point vortex 
of the order parameter is introduced as 



Ai +iA 2 = |A(r 



Jnrj) 



(2) 



where n is the vorticity, and (f> is the azimuthal angle, 
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and |A(r)| describes the radial variation of the ampli- 
tude. It has been shown via a direct calculation^ - that 
there are exactly n number of zero modes for vorticity 
n. Subsequently, the correspondence between the vortic- 
ity and the number of the zero energy states, has been 
expressed as a Z index theorem by Weinberg^ 2 -. This 
is an index defined on the open spaced like Callias in- 
dex theorem^, unlike the celebrated Atiyah-Singcr index 
theorem 2 ^, which is defined on the compact manifold. 
There is an additional unitary matrix T 5 = r 1 r 2 r 3 r4, 
which anti-commutes with Hjr, and ensures its spec- 
tral symmetry. All together the five mutually anti- 
commuting T-matrices close the Clifford algebra of the 
four dimensional matrices^. As a consequence of the 
spectral symmetry, the zero energy states become eigen- 
states of r$, with eigenvalues ±1. Therefore T$ defines 
the chirality of the zero modes. The Z index theorem 
states that 

N+ - A_ = n, (3) 

where N± are respectively the number of zero energy 
states with chirality ±1. In the Altland-Zirnbauer 
classification scheme, this Hamiltonian belongs to the 
class BD I 27 ' 28 . The presence of the gauge field is 
not important for the number of zero modes or the 
associated index theorem. However, the existence of 
the gauge field is of paramount importance to realize 
deconfined vortices. 

The Jackiw- Rossi Hamiltonian Hjr [a] can describe 
both insulating and superconducting states. For an in- 
sulator, aj's are the components of the chiral gauge field. 
In contrast, for a superconductor a^s describe the elec- 
tromagnetic vector potential. For the insulating states 
of the charged fermions, the electromagnetic U(l) sym- 
metry is preserved and the Hjr [a] can be augmented by 
the electromagnetic vector potential A, without altering 
the BDI chiral symmetry^ 9 .. The modified Hamiltonian 
is 

H JR [a,A] = ]T r j {-id j +iT 3 T 4 a j +A j )+A 1 T 3 +A 2 T A , 

.(4) 

The vector potential A only changes the wave functions 
of the zero modes, without affecting its numbe r 29 ' 30 . In 
particular for the insulators, when the vorticity is one 
(n = 1), there exists single state at zero energy. This 
leads to the charge fractionalization, which is relevant 
for the spinless fermions on the honeycomb lattice with 
an underlying vortex of the Kekule bond density wave 
orde r 12 ' 31 " —. The axial magnetic field (time-reversal 
preserving) can be introduced on the honeycomb lattice 
by deliberate wrinkling of the flako^i. For graphene, 
due to the spinful nature of the fermions, the charge 
fractionalization is lost. Instead, there is a competing 
antifcrromagnctic order in the Kekule vortex core^ 
The fcrmion zero modes also occur for the vortices of 
the easy plane antiferromagnetic and quantum spin Hall 
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The fermion zero mode for the s-wave pairing of the 
topological insulator's surface states is also described 
by Hjr [a] (by ignoring the chemical potential and 
the Zeeman coupling). Due to the BdG nature of the 
quasi-particles, the single zero energy mode corresponds 
to a Majorana fermio n 37 ^ 38 . The superconductivity 
on the helical surface states of a topological insulator 
has been realized by the proximity effect^ 9 ., and the 
possibility of realizing the Majorana zero mode is also 
being actively pursued^ .. Similar considerations can be 
applied for the superconducting states on the honeycomb 
lattice 4 ^—. However, the minimal representation of 
spinful Dirac fermions in graphene is 8-dimcnsional and 
to accommodate the pairings, the Nambu's doubled 
representation becomes 16-dimensional. In this context 
the associated Hamiltonian is a direct sum of multiple 
copies of the four dimensional Jackiw-Rossi Hamilto- 
nian, shown in Eq. ([T]). The existence of multiple zero 
energy modes even with n = 1, ruins the possibility 
of realizing Majorana fermions 4 ^. But, the presence of 
the zero modes leads to an interesting interplay of vari- 
ous competing order parameters inside the vortex core^ 

So far, we have been ignoring the effects of the finite 
chemical potential and the Zeeman coupling on the zero 
modes in a superconductor. Due to the inevitability of 
these couplings in a real system, it becomes important 
to ask which generalization of H JR [a] still supports zero 
modes. The answer to this question was provided in a 
recent work by Herbut and Lu 4 ^. It has been argued 
that one can introduce two additional parameters to the 
original Jackiw-Rossi Hamiltonian in Eq. ((T|), which to- 
gether can support zero energy mode. The generalized 
Jackiw-Rossi Hamiltonian then takes the form 

H'll = Hjr [a] + ir 3 r 4 A + iY 1 Y 2X . (5) 

It has been shown in Ref. [46| that there exists an anti- 
unitary operator (A) , which anti-commutes with the total 
Hamiltonian Hg^, guaranteeing its spectral symmetry. 
In the same work authors have showed that a point vor- 
tex of vorticity ±1 can bind one zero energy state only 
if Aq + A 2 > x 2 , where A = \A(r oo)|. The physi- 
cal meaning of the parameters (A, x) are again represen- 
tation dependent. For example, if we consider spinless 
fermions in monolayer graphene, A corresponds to a chi- 
ral chemical potential, whereas x to the Haldanc anoma- 
lous mass4£ On the other hand, for s-wave paring of the 
surface states of topological insulator, they respectively 
correspond to the ordinary chemical potential and the 
Zeeman coupling. It has been argued recently that the 
generalized Jackiw-Rossi Hamiltonian can support a sin- 
gle zero energy mode only if the vorticity is od d 48 ' , and 
all the states are at finite energies for the even vortic- 
ity (except at A = ±x)^°"- Since, the spectral symmetry 
with respect to T§ (unitary) is now lost, the concept of 
chirality and Z index becomes moot. Therefore, the gen- 
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cralized Jackiw-Rossi Hamiltonian has a Z2 index for the 
generic parameters (A 7^ ±x): which can be stated as 

M=\ (l-(-l)") e(Ag + A 2 - x 2 ), (6) 

where M and n are respectively the number of zero 
modes and the vorticity. A similar Z2 index also emerges 
for the two dimensional p + ip supcrconductor a 11 ' 51 , 
which has been argued to be a non-relativistic limit of 

ttip HJR 50,52 
gen 

In this paper we will mainly focus on the time rever- 
sal symmetric, fully gapped, momentum independent 
superconducting states of three dimensional Dirac 
fermions, which are realized in many narrow gap and 
gaplcss semiconductors. In particular, there has been 
a surge of interest in the superconducting states of 
e.g. copper intercalated bismuth selenide Cu^B^-Se^, 
indium doped tin tclluride Sni- x In x Te£^ Apart from 
the regular s-wave pairing, the three dimensional Dirac 
nature of the quasi-particle also allows the possibility 
of a fully gapped topologically non-trivial odd-parity 
pairing^ Even though the pairing symmetry in these 
materials has not been established^ it is known that 
the paired state is fully gapped and of typc-II nature. 
This motivates us to study the possible dispersive 
modes along the vortex core in the mixed phase of these 
materials, and for concreteness we restrict ourselves in 
the dilute vortex limit (close to H c i). 

Now we provide a synopsis of our main findings in this 
paper. 

1. For the -ff/ e „, we explicitly demonstrate the 
algebraic reason for the absence of the zero modes 
in the case of even vorticity. In addition, we 
reconfirm the existence of the single zero mode for 
odd vorticity. 

2. We first consider the pairing of three dimensional 
massless Dirac fermions in the absence of the 
chemical potential and all other fermion bilinears 
in the particle-hole channel. When projected to 
the x — y plane (k z — 0), the Hamiltonians for both 
s-wave and the topological odd parity pairings 
map to two copies of Hjfj[0}. Therefore, the Z 
index theorem of Hjr[0] governs the number of 
Majorana zero modes. 

3. When we incorporate the chemical potential, the 
Hamiltonian for both pairings map onto two copies 
of Hg^. Consequently, a pair of Majorana zero 
modes are only found for odd vorticity. In addition, 
if we introduce the Dirac mass and the Zeeman 
couplings, the odd parity topological pairing still 
leads to two copies of H J g f n . On the other hand, 
the effective Hamiltonian for the s-wave pairing 



docs not map onto i?g e „ in the presence of any 
chiral symmetyry breaking perturbations in the 
particle-hole channel (for example scalar Dirac 
mass). Consequently, the zero modes are absent 
for s-wave pairing under generic perturbations. In 
contrast the zero modes for the the topological 
superconductors are more robust. 



4. For the topological superconductor, k z causes mix- 
ing between the two Majorana modes and converts 
these into a dispersing complex fermion. We pro- 
vide a symmetry argument to explain why k z does 
not cause any mixing between the zero modes and 
the finite energy Caroli-de Gennes-Matricon states. 
Due to this reason, the degenerate perturbation 
theory for k z in the zero energy subspace provides 
an exact answer for the dispersion relation. We 
also demonstrate this through some exact solutions. 



5. In the concluding sections, we consider the 
generalization of the zero mode problem in the 
context of axial superfluid. We show that the 
fermion zero modes can be found in the presence 
of the axial chemical potential and the third 
component of the axial vector, for odd vorticity. 
This leads to an interesting modification of the 
Callan-Harvey mechanism of anomaly cancellation. 



The rest of the paper is organized as follows. In the 
Sec. II we demonstrate the Z2 index associated with the 
generalized Jackiw-Rossi Hamiltonian. In Sec. Ill, we 
explicitly present the mapping of the three dimensional 
Dirac Hamiltonian with pairings to two copies of the 
generalized Jackiw-Rossi Hamiltonian under appropri- 
ate conditions. Sec. IV is devoted to establish the 
invariance of zero energy sub-space under the influence 
of the momentum along the vortex line. In Sec. V, 
we exemplify this claim by computing the dispersive 
modes with energies E = ±k z , as well as the vortex 
zero energy modes for particular choice of parameters. 
There we show that the dispersive modes can indeed be 
constructed from two dimensional zero energy modes 
by performing a perturbation calculation in terms 
of k z . In Sec. VI, we summarize our main findings 
regarding the Z2 index of HjfJ* and its consequences 
for superconducting states of three dimensional Dirac 
fermions. In addition we propose a generalization of the 
Callan-Harvey model of axial superfluid, and discuss 
the consequences Z2 index on anomaly equations. We 
relegate the details of the zero mode calculations at the 
special points A = ±% to the Appendix A. 
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II. Z 2 INDEX OF GENERALIZED 
JACKIW-ROSSI HAMILTONIAN 

Upon setting set A = x = in Eq. (0, the resulting 
Hamiltonian conforms to the one studied originally by 
Jackiw and Rossi for point vortices^ It was shown in 
Ref. and Ref. that for arbitrary vorticity (n) of 
the point vortex, there exists precisely n number of zero 
energy states. This problem belongs to class BDI in the 
Altland-Zirnbauer classification . 27 ' 28 The generalized 
Jackiw- Rossi Hamiltonian in Eq. ([5]), may as well 
support states at zero energy depending on the vorticity 
and the relative strength of the parameters, A, x an d the 
asymptotic value of the order parameter as r — > oo (A ). 
Possible existence of the zero energy mode is ensured 
by the spectral symmetry of the Hamiltonian, H^. 
Such symmetry of Hg^ is generated by an anti-unitary 
operator, namely A = UK, where U is the unitary 
operator and K is complex conjugation . 46 ' 57 Without 
any loss of generality, one can commit to a representation 
where Ti,r2 are real and T3,r4 are imaginary 26 ' 46 ' 58 
In that representation, U is the identity operator, 
and the anti-linear operator A is simply the complex 
conjugation. Furthermore, if there exists any state at 
precise zero energy in the spectrum of Hg^, it needs to 
be an cigenstatc of A with eigenvalue ±1. 



A. Spectral symmetry and zero modes 

Let us now focus on the zero energy modes of the 
Hamiltonian H^^, in Eq. Since all the four di- 

mensional representations of mutually anti-commuting 
matrices are unitarily equivalent, for our convenience we 
choose to work with 5 ^ 

Ti = -0-3 (8 <7i, r 2 = Co <s> 0-2 , r 3 = el (g) <7i, 
r 4 = cr 2 <g> (Ji , r 5 = rir 2 r 3 r 4 = a ® o- 3 . (7) 

Let us define a four component spinor as 

V> T (£) = {ui,vi,u 2 ,v 2 ) (x), (8) 
and here we wish to solve 

Kn <K*) = - ( 9 ) 

In this representation, the anti-linear operator, which 
anti-commutes with Hg^ and ensures its the spectral 
symmetry is 

A = iT 2 T 3 K = (o-iig) cr 3 ) K. (10) 

This anti-unitary operator leaves the zero energy sub- 
space invariant. Hence, the zero energy state needs to be 
an eigenstate of the operator A, with eigenvalue either 
+1 or —1, implying 




Upon imposing the constraint on the spinor components 
with the + sign in last equation, the four coupled differ- 
ential equations for the zero energy mode reduce to only 
two, which read as 

(-i)e i *(dr + ldAvZ + A r e- m *v 2 + (\ + xK 2 = 

{i)e-^{d r ~ l -d^ju* 2 + A r e- m4 u 2 -{\-x)v* 2 = 0. 

(12) 

For the ease of calculation we here chose to work with 
an underlying anti- vortex. Our discussion is equally ap- 
plicable for point vortex. The above set of equations has 
been derived for a = 0. However, the structure of these 
equations are not qualitatively affected by a, and we will 
explicitly account for the gauge field in the subsequent 
sections. From the above equations, the imaginary factor 
can be removed by redefining the spinor components as 

«*(£) = e**V*{Z), u* 2 (x) = e~^U*(x). (13) 

Next we find the zero mode solution of the generalized 
Jackiw- Rossi Hamiltonian, shown in Eq. ([9]), with 
underlying point vortex of even and odd vorticitics. One 
may have taken the zero energy state, as an eigenstate 
of A with eigenvalue — 1. However, this corresponds 
to a phase rotation by exp(i^r 5 ), which is not an 
observable^ 



B. Even vorticity: n — 2s 

Let us first consider even-vortex, i.e. n = 2s, where s is 
a positive integer. With even vorticity single-valued so- 
lution of the zero modes can only be found assumin g 41 ' 60 

V(x) = e^A(r)+e"^ 5l (r); 

U(x) = e**/ 2 (r)+e <9 *<fe(r). (14) 

The consistent solutions of Eq. ((9]), with the above ansatz 
can be obtained upon imposing the following constraints 
over the angular momenta 

l + m = 2s + l;p + q = 2s-l;p = l-l;q=2s-l. (15) 

One should notice that the above constraints among dif- 
ferent angular momenta channels of V and U arise only 
when A and/or x are nonzero. Otherwise, U and V are 
decoupled from each other, same as in the original work 
by Jackiw and Rossi^ - Upon imposing these constrains 
over the angular momenta, we obtain the following cou- 
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d r 4 

d r - 

d r - 

d r 4 

where 



pled differential equations for the zero energy modes 

/ 1 (r)+A r5l (r) + (A + x)/ 2 (r) = 0, (16) 

/ a (r) + A r92 (r) - (A - X )/i(r) = 0, (17) 

g 2 (r) + A r / 2 (r) - (A - ^(r) = 0, (18) 



r 

I - 1 
r 

2s - I 



2s + 1-1 



g 1 (r) + A r f 1 (r) + (X + X )92(r) = 0, 

(19) 



1 = 0, ,(2s-l). 



(20) 



Therefore, with even vorticity, n = 2s, there may 
exist n number of zero energy states. It is quite 
challenging to obtain the full analytical solution of 
above four equations. However, we here show that the 
existence of the zero energy modes can be established 
by studying the asymptotic behaviors of the solutions: 
1) Near the origin, r — > 0, where the pairing order 
parameter A r vanishes smoothly, and 2) far away 
from the origin, r — > oo, where the A r approaches a 
constant value Ao and all the terms proportional to 1/r 
in Eq. (|T5|) - (fT5|) can be neglected. Next, we proceed 
to hnd the asymptotic solutions of the zero energy states. 

Near origin (r — > 0): Neglecting the contribution from 
A r , we arrive at the two sets of coupled differential equa- 
tions. One of them reads as 

^■ + ;)/i < (r) + (A + x)/ 2 < (r) = 0, 



r — 



I - 1 



f?(r)-(\-x)fHr) = 0. (21) 



Here, we denote all the radial functions in the vicinity of 



the origin as A" 



where X = f 1 , f 2 , gi, 92- The other 



set of the coupled differential equations is 
2s + 1 - V 



Or 



r 

2s 



9?{r) + {\ + xm{r) 



5 2 < M-(A-xK < W 



0. 

0. (22) 



Solutions of the first set of equations (Eq. 

r<f^ -I °l J VTU> V gZZ3 if A > x ' 



arc 



/r(r) 



dl^r^yj^] ifA< X , 



(23) 



where C; is an arbitrary constant. Jk and Ik are re- 
spectively the Bessel and the modified Bcssel functions 
of first kind of order k. From the solution of fi{r), one 
can immediately find (r) from 



fHr) = 



1 



A + X 



9 r + ~). 



The second set of coupled equations (Eq. (|22l) ) yields 



(25) 



Here C; is also an the arbitrary constant. And the re- 
maining function, g^ (r) can be found using 



1 



A + x 



dr 



2s + l-l 



g<(r). (26) 



Therefore all the four radial functions in the vicinity 
of the origin is defined in terms of one (1) arbitrary 
constant. 

Far from origin (r — > 00): Next we proceed to find 
the radial dependence of these four functions far from 
the origin. Neglecting all the terms proportional to 1/r 
in Eqs. (fHfl) we obtain a new set of four coupled 

differential equations 

a r />(r) + Ao.^ (r) + (A + x)/ 2 > M = 0, 

d r g>(r) + A />(r) + (A + x )g>(r) = 0, 

d r fi(r) + A g>(r) - (A - x )/>(r) = 0, 

d r g$(r) + A />(r) - (A - x )g?( r ) = 0. (27) 

Radial dependence of all the functions far away from the 
origin is denoted by X y (r), where X = /1, / 2 , 51, g 2 . No- 
tice that far away from the origin the differential equa- 
tions are independent of the angular momenta (I, m, p, q). 
The above four equations reduce to two sets coupled dif- 
ferential equations in terms of new variables, defined as^i 

F±(r) = />(r) ± g>(r); and G ± (r) = />(r) ± <? 2 >(r). 

, (28) 

In terms of these variables the set of equations in Eq. 
(f2"T|) away from the origin becomes 

d r F ± (r)±A F ± (r) + (X + X )G ± (r) = 0, (29) 
d r G ± (r)±A F ± (r)-(\- X )F±(r) = 0. (30) 

The solution of these new equations can in general be 
written as 



F + =^2,C a exp (a£r) , G+ 

<T=± 

F_ = Y,Ca exp (c£r) , G_ 



^^exp^r), 
C* exp (a^r) , 



<T=± 



(31) 



where 

of = -A + cr^fx^y?; a G a = A + a^ X T ^, (32) 

with a = ±. The arbitrary coefficients appearing in the 
solutions are related according to 



C' 



V X + A 



'X_A 
X + A' 



(33) 
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We are interested only in those solutions which decay 
exponentially as r — > oo, so that they are normalizable. 
Depending on the relative strength of Ao, A and x, the 
solutions can take different forms. We study each of 
these cases separately below. 

If > x 2 - A 2 < Ag, a*± < 0, but al > 0. Therefore, 
normalizibility of the solutions demands _F_ = G_ = 0. 
In terms of the original functions the solutions are 

/f» = 9Ur) = \ (c+ cxp (o£r) + C_ exp (a F r) 



f?(r)= 9 >(r)= 1 -J X ^ ' r 



2V X + A 



C+ exp (a 1 7 



-C_ cxp (a F r) ^ . (34) 



If on the other hand, x 2 — A 2 > A 2 ,, a_ G < 0, whereas 

a+ G > 0. Hence, for normalizable solutions C + = C + = 
0, and we obtain 



1 



C_ exp (o^r) + CL exp (oPi 



9?(r) = -( C_exp(c/» -C_exp(a G r) 



(35) 
(36) 



gUr) 



2VX + A 



C_ exp (afr) — C- exp (a G r) I , 



(37) 



1 / X -A 

2 V X + A 



C_ exp (a_r) + C_ exp (a_r) 



(38) 



Therefore, irrespective of the mutual strength of A, x an d 
Ao, the radial part of the zero mode solutions far away 
from the origin is always defined by two (2) arbitrary 
constants. 

Boundary conditions: To obtain self consistent solu- 
tions of the zero energy states, the spinor components 
need to satisfy the boundary conditions at a particular 
point r = £(say). Since the solutions in the asymptotic 
regions arc obtained by solving a second order differ- 
ential equations, for each functions, we need to match 
values and the first derivatives of X < {r) and X > (r) at 
r = £, where X = fi,gi, f2>92- However, for example, 
when we impose these boundary conditions over /2 (r) , it 
remove the arbitrariness from two out of the three con- 
stants, cither (C + ,C_,C;) or \ C-,C-,CiJ, depending 

on whether x 2 — A 2 < or x 2 — A 2 > Aq , respectively. In 
either situation, such elimination immediately removes 
the arbitrariness of the constants in the solution of (r) 
as well. Therefore, with two fixed constants in g^(r) 
and one arbitrary constant in the definition of gf(r) 1 it 
is impossible to satisfy two of its boundary conditions. 



Similarly, upon imposing the boundary conditions over 

/i (r) , two out of cither (C+ , C*_ , Q ) or {(J- ,C-,C^ get 

fixed, depending on x 2 - A 2 < or x 2 — A 2 > A§, re- 
spectively. In terms of those fixed constants it is once 
again impossible to satisfy both the boundary conditions 
for 172(f)- Hence, one can conclude that for generic pa- 
rameters (A 7^ ±x)j there exists no zero energy mode for 
Hgln' when the vorticity is even. 



C. Odd vorticity: n = 2s + 1 

Next we consider point vortex with odd vorticity, and 
take n = 2s + 1, where s > is an integer. The single- 
valucdncss and the consistency of the zero mode solutions 
immediately imply that out of 2s + 1 choices for the zero 
modes, there is precisely one solution for which we can 
choose 

V(S)=e a +f(r); U{x) = e»*+g(r), (39) 

and the angular momenta I and p are related as 

l = s + l=p+l. (40) 

With these constraints over the angular momenta, the 
coupled differential equations for f(r) and g(r) read as 

d r + ^jf(r)+A r f(r) + (X + X )g(r) = 0, 

d r + ^jg(r) + A r g(r)-(X-x)f(r) - 0.(41) 

Even though, we can exactly solve these two coupled dif- 
ferential equations, it is worth analyzing the existence of 
the zero mode solutions with the above ansatz, from their 
asymptotic behaviors. Later we will present the complete 
solution and show that these two approaches complement 
each other. In the vicinity of the origin (r — > 0), upon 
dropping the contribution from A r , we obtain 



Ql^r^x^V} ifA<X, 



and g < (r) can be found from the relation 

g < (r) = -(A-r) (dr+ l -)f<(r) 



A + x 



(42) 



(43) 



At large distances (r — > 00), upon setting A r = 
Ao (constant) and after dropping the terms proportional 
to 1/r, we obtain 

d r f > (r) + A f > (r) + (X + x)g > (r) = 0, 
d r g > (r) + A g>(r)-(X- X )f > (r) = 0, (44) 

which can be solved by using the ansatz f > (r) — 
c exp (ar) , where 



-A ± v 7 ? 



A 2 



a±. 



(45) 
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However, only for x 2 — A 2 < and consequently A 2 +A 2 > 
% 2 , the solutions f > (r) and g > (r) are expressed in terms 
of two arbitrary constants. Then the solutions at large 
distances are 



at large distances the modified Bessel functions grow ex- 
ponentially 



'ar) oc 



(53) 



9 > (r) 



ix-x 

A + x 



^2 vC a Q a (A ,x,\r), 



<7=± 



/>(r) = VCA(A„, X ,V), 



<T = ± 



where 

<2<x(A ,x, A,r) 



exp -A r 



^x/A 2 



(46) 



(47) 



Therefore, upon imposing the boundary conditions, we 
mentioned in the previous sub-section (say r = £), two 
out of three arbitrary constants for each function get 
fixed, while the remaining one is then set by the nor- 
malization condition. If the parameters are such that 
Aq + A 2 < x 2 j there is only one exponentially decay- 
ing solution, cx cxp(a_r). Therefore, each function is 
defined in terms of one arbitrary constant at large dis- 
tances. Therefore, we have two arbitrary constants and 
two matching conditions to satisfy. After satisfying the 
matching conditions, there is no more arbitrary constant 
is left to set the overall normalization. Hence, we can 
have atleast one normalizable zero mode of H J g f n for odd 
vorticity and 



Therefore, normalizable zero energy modes can only be 
found when, A 2 , + A 2 > x 2 , identical to the one we found 
in Eq. (|48[) by analyzing the asymptotic solutions. On 
the other hand, when A > x, f(r) is defined in terms of 
the Bessel functions of the first kind, as shown in Eq. 
(|5Tj) . The condition in Eq. (|4"5)l is then trivially satisfied 
and we always find a normalizable zero mode. 

Besides the above zero mode solution, there are addi- 
tional 2s possible ansatz similar to that in Eq. (fT4|) . The 
angular momenta satisfy the following constraints 



l+m = 2s+2,p + q = 2s, p = l-l,q= (2s + 1) 



(54) 



A 2 + A 2 > x 2 . 



(48) 



Upon imposing this set of constraints over the the angular 
momenta, one set of coupled differential equations for the 
functions fi(r) and f£~{r), in the vicinity of the origin 
assumes the identical form as in Eq. (f2"Tj) . Their solution 
can readily be found from Eq. (|23|) . (|24|) . The other set 
of equations in terms of g^ (r) and g^~ (r) is 



9r + 25 + r 2 l \f (r) + (A + x)g 2 < (r) = 0, 

d r + l —^±)g<(r) - (A - X )g<(r) = 0, (55) 



Next we show that this condition can also be obtained 
from the exact solution of the zero energy mode. 

The exact solution of the zero energy states from Eq. 
(gT]) can be found upon assuming 



yielding 



?Mf3 i [ A ! X ' (56) 
I °i I V2l+2^i[ r VX 2 - A 2 J if A < x- 



f(r) = f(r)e~f° Ar ' dr ', g(r) — g(r)e~ & <V*-', (49) One can then find gf (r) from the relation 



The functions f(r) and g(r) then satisfy the following 
equations 



A + X 



Or 



2s + 2-l 



(57) 



Or 



d r — 



f(r) + (A + X )g(r) = 0, 
ff(r) + (A + x)/(r) = 0, 



yielding 
and 



dljyr^^] ifA<x, 



A + x 



(50) 



(51) 



(52) 



If wc assume x > A, then the solutions are expressed in 
terms of modified Bessel functions of first kind. However 



Hence, at small distances there is one arbitrary constant 
for each of /i(r), f2(r), gi{r) and g^ (r). On the other 
hand, we have shown in the previous section that the 
large r behavior is independent of the angular momenta. 
Hence, far away from the origin the radial dependence 
is captured by the Eq. ([Mf if x 2 - A 2 < or Eqs. 
(|35 ]) -(|38 l) if x 2 - A 2 > A§. Therefore, following the 
discussion on the matching conditions at r = £, in 
the previous subsection, we can argue that there is no 
zero energy mode with multiple angular momenta ansatz. 



D. Effects of gauge potential 

In the above derivation, we have neglected the orbital 
effects of the gauge potential. It can be introduced, for 
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example, considering a simple profile of the magnetic 
field. Let us assume, that the magnetic field (applied 
in the z-direction, perpendicular to the plane of the vor- 
tex) is finite and constant only within a distance r < £, 
and vanishes for r > ^i 46 ' 62 Then in the symmetric gauge, 
one can choose = t£z when r < £, and = l/2r for 
r > £. With such a profile of the gauge potential, the 
ultraviolet and the infrared asymptotic behaviors of all 
the functions (/i, /2, gi, 52, /, <?) remain unchanged. Only 
significant effect will be at intermediate distance r ~ £. 
Hence, our derivation for the Z2 index for the general- 
ized Jackiw- Rossi Hamiltonian remains valid, even in the 
presence of the gauge field. To further demonstrate this 
assertion, we explicitly solve the vortex zero mode with 
this profile of the gauge field in Sec. V. 



describing a point vortex defect by unitarility rotating 
Hjr in Eq. ([I} by iT^T^, which changes the relative 
sign between Ai and A2. When we perform the same 
operation on the generalized Jackiw-Rossi Hamiltonian 
in Eq. (J5)) it changes the sign of the U(l) gauge field 
a — > —a, and takes A — > —A. Hence, n number of zero 
energy modes with underlying point vortex appears 
in the spectrum when A = +x, whereas the Z2 index 
remains unaltered for A = — X- With an underlying 
vortex when A = +x, all the n number of zero modes are 
eigenstate of T5, however, with eigenvalue —1. For the 
detail solutions at this special values of two parameters 
A and x, readers are referred to the Appendix A. 



E. Statement of Z2 index 



III. THREE DIMENSIONAL MASSIVE DIRAC 
HAMILTONIAN WITH GAPPED PAIRINGS 



After going through the above arguments, we can 
formally present the statement of the Z2 index associated 
with the zero modes for the generalized Jackiw-Rossi 
Hamiltonian. For generic values of A ^ ±x, which 
satisfy A 2 + A 2 > \ 2 , there exists one normalizable 
zero mode when the vorticity is odd, while with even 
vorticity all the states are placed at finite energies, in 
accordance with Eq. ©. On the other hand, there 
are no normalizable zero modes for any vorticity, if 
Aq + A 2 < x 2 ■ We notice that the normalizability 
at infinity is governed by the uniform H 9 ^ . There 
is a topological phase transition or band inversion at 
Aq + A 2 = x 2 , and this is the reason for the absence of 
the normalizable zero energy state on one side. 



Next we focus on momentum independent, time rever- 
sal symmetric, gapped paired states of the three dimen- 
sional Dirac fcrmions. In three spatial dimensions the 
Dirac quasi-particles can pair into two fully gapped su- 
perconducting states. One of them is trivial s-wave pair- 
ing, whereas the other one is parity-odd and topologically 
non-trivial. To study the nature of these paired stated in 
the mixed/ vortex phase, let us define an eight component 

Nambu-Dirac spinor as ^ - 
*J(+k) = * T (k) and *T(- 



^J(+k),^J(-k) , where 
k) = fy p (k), otherwise 



4(fc),c+(fc),cf (k),cj(k) . (58) 



However, there is an exception to the existence of 
the Z2 index at A = ±x^ For these special values, 
the condition A 2 , + A 2 > x 2 is trivially satisfied, which 
guarantees the normalizability of the zero modes. In 
the above calculation, we have assumed an anti-vortex 
configuration of the order parameter (A). With an 
underlying anti-vortex it can be shown that there exist 
precisely n number of zero energy states for arbitrary 
n, when A = — x- Moreover, all the zero energy states 
are eigenstates of the the chirality operator r$, with 
eigenvalue +1. On the other hand, for A = +x, once 
again we recover the Z2 index for the generalized Jackiw- 
Rossi Hamiltonian. One can achieve the Hamiltonian 



cf corresponds to the annihilation operators for the even 
and the odd parity states, respectively, with the spin pro- 
jections s =f, 4- The three dimensional massive Dirac 
Hamiltonian in the presence of trivial s-wave (A s ) and 
odd-parity topological (At) parings takes the form 



H D = ^2^(k)H gen 

k 



[A S) A T }*(k). 



(59) 



In order to preserve the time reversal symmetry we do 
not choose any relative phase between A s and At. In 
the announced 8-dimcnsional Nambu-Dirac basis (VP), 



H gen [A s , At] takes the form 



Hg en [A s , At) = (r ® ai)(k x - r 3 ® I 4 a x ) + (r 3 ® a 2 )(k y - r 3 ® I A a y ) + (to (8 a 3 )(k z - t 3 ® I&a z ) + (r 3 ® h)fi 
+ (t 3 ® P)m k - (t 3 <g) iaia 2 )h + - (t 3 ® ia^T)h- - (t 2 ® a 2 )A^ + (t\ ® a 2 )A^ 
- (r 2 ® iaia 3 )Af + (n ® iaia 3 )Af , (60) 

I 

where a is the electromagnetic gauge potential, and the with x = T,s, and n is an integer. Then n counts the vor- 
complcx pairing order parameters are defined as ticity and the Hamiltonians H gen [A s , 0] and H gen [0, At] 

A x = (A* A£) = |A K |(cosn0,sinn0), (61) 
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respectively correspond to line-vortex defect with un- 
derlying s-wave and topologically non-trivial odd parity 
pairing. h± are respectively the symmetric and the anti- 
symmetric combinations of the Zeeman couplings of the 
even (hi) and the odd (h 2 ) parity bands, namely 



h± = ^\hi±h 2 \ 



(62) 



fi is the chemical potential and nik is the Dirac mass. 
In what follows we set = m = constant. The four 
dimensional Hcrmitian matrices are defined as 











0-2 ] 




( 


0-3 \ 


V CT i 




I tx 2 


o J 


, "3 = 1 


V 0- 3 


) 



a Q 
-<T 



— ia 
i<7Q 



(63) 



which together complete the Clifford algebra of five 
mutually anti-commuting four dimensional matrices. 
Here, o~o is the two dimensional unit matrix and 3 arc 
the Pauli matrices. The other set of two dimensional 
matrices, {tq,t}, operate on the Nambu's index. Here 
we have ignored the anisotropy in the Fermi velocity 
arising from the underlying crystallographic structure, 
and set v x = v y = v z = v = 1j22 Next we cast the pairing 
Hamiltonians in the k z = plane with underlying point 
vortex defects of the topologically non-trivial odd parity 
pairing as orthogonal sum of two copies of the general- 
ized Jackiw-Rossi Hamiltonian, under generic situation. 
Such mapping is shown to be true for the s-wave pairing, 
however only if there is no chiral symmetry breaking 
perturbations, e.g., to, 



and 



H 



1 = I5U1K - /3"f 5 a 2 k v - i/3a 2 A% - 
+ ^"375 (m - h + ) + /3 (/x - /i_) 



a 2 Al 



(67) 



The matrix multiplying k z is Ai^ = a 2 <E> (10(1013). Here, 
we have defined as new matrix 



75 = 



0-0 
<J 



(68) 



which anti-commutcs with /3 and T, but commutes 
with a\,a 2 ,ct3. Both H^ L ^ dR in Eq. (|66p . and Eq. 
(|6"T)) assume the form of the generalized Jackiw-Rossi 
Hamiltonian, shown in Eq. We need the following 

identification of the parameters \ = m + A = // + /i_ 
and x = m — A = // • 



for H% L 



h_ for H% R . 



Therefore, with a point vortex defect of underlying 
odd-parity topological paring one ends up with two Ma- 
jorana modes, in the presence of generic perturbations. 
However, these Majorana modes can only be found when 
the vorticity of the point vortex is odd, as we have shown 
in the previous section. Since the magnetic field gets 
screened beyond the core of the vortex, the Zeeman cou- 
plings h± are finite within the vortex core. Hence, two 
normalizable Majorana modes can be achieved only when 



2 ^ 2 

(J, > 771 . 



(69) 



s-wave pairing 



A. Odd-parity topological pairing 

To perform the above mentioned exercise, it is worth 
redefining the 8-componcnt Nambu-Dirac spinor as 



c+ c7, (c+) f , (c7) f , c+ cf , (c+) f , (cf ) f ] (x). (64) 



In 
H 



this new basis the part of the Hamiltonian, 
gen [0, Ay] describing the point vortex in the rry-planc 



(i.e. with k z = 0), is completely block-diagonal, whereas 
the k z part is block off-diagonal. For simplicity, let us set 
all the orbital components of the gauge potential to zero. 
The total Hamiltonian with only the topological paring 
then takes the form 



H gen [0,A T ]^H^ r = 



H'l 



H 



dR 



+ k z M T z , (65) 



where 



ttuL 

£lj> 



j 5 aik x + (3j 5 a 2 k y + i(3a 2 A R + a 2 A^ 
pa 3 j 5 (m + h + )+0((x + h-), (66) 



A similar exercise can also be performed with an un- 
derlying s-wave pairing, which also assumes the form of 
the generalized Jackiw-Rossi Hamiltonian, when m = 
and h- = 0. Let us now rewrite H gen [A s , 0] in the basis 



c+ cr, (c+)t, (cf )t, c+ cf , (c+)t, ( c -)t] (£). (70) 



The eight dimensional Hamiltonian, H gen [A s , 0] then 
takes the form 



H gen [A a , 0] -> fl?" = (^ S " L © iff M + fc z A^, (71) 

similar to the Eq. (|65|) . The diagonal blocks of Hg° r are 

Hf L = 7s (ai^ + a 2 k y ) + Af a 3 + Afia 3 /3 + M/3 

+ 75^3^+ — */3aia 2 TO + / 4 /i_, (72) 

and 

fff* = 7 5 (aifctt-a 2 A; ! ,)-Afa3-A^a3/3 + M/3 

— 7503/1+ — z/3«ia2TO — I^h—. (73) 

Here as well the k z appears as block off-diagonal clement, 
and the matrix multiplying k z is M. s z = a 2 ® (iaia.3), 
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identical to Ai^ . Therefore, in the absence of the 
Dirac mass (m) and the anti-symmetric combination 
of the Zeeman coupling (h-), both H^ L and Hf R are 
equivalent to the generalized Jackiw-Rossi Hamiltonian, 
shown in Eq. ([5]), where A = /x, x = h+ f° r both H g eJ dR . 
Hence, in the absence of any chiral symmetry breaking 
perturbations, the point vortex of underlying s-wave or- 
der can also support two Majorana zero modes, when the 
vorticity is odd. The zero energy modes in the absence 
of /i, h + , are also the eigenstates of /3j 5 = — er 2 ® cr , 
with definite chirality. This matrix defines the chirality 
of massless three dimensional Dirac fermions, and 
should not be confused with T§ — <tq ® (T3 — ia.2Ct\, 
which defines the chirality of the original Jackiw-Rossi 
Hamiltonian in Eq. (fTJ). Any chiral symmetry breaking 
perturbation of the three dimensional Dirac fermions, 
such as Dirac mass (to) and h- cause mixing among 
these two states and the spectrum becomes gapped. 



IV. PERTURBATION THEORY FOR LINE 
VORTEX ABOUT k z 

In this section we show that the Z 2 index for the 
fermionic zero modes bound to the point vortex, also 
dictates the number of one dimensional dispersive modes 
along the line vortex in three spatial dimensions. The 
momentum along the vortex core (k z ) is shown to leave 
the zero energy sub-space of the underlying two dimen- 
sional Hamiltonian (multiple copies of Hg^) invariant. 
Consequently, a first order perturbation theory for k z is 
exact and the dispersive modes along the vortex core are 
the linear combinations of the vortex zero modes. We 
here prove this statement with underlying topologically 
non-trivial odd-parity and s-wave pairing separately 



A. Topological pairing 

Let us first consider the generalized Dirac Hamilto- 
nian with the topological pairing, H gen [0, At]- Wc 
have shown in the previous section that H gen [0, At] 
is unitarily equivalent to two copies of the generalized 
Jackiw-Rossi Hamiltonian, when k z = 0. Therefore, 
H gen [0, At; k z = 0] hosts Mojaorana zero modes for odd 
vorticity, under generic situation. These two zero en- 
ergy states constitute a two dimensional basis, which re- 
mains invariant by any operator that commutes or anti- 
commutcs with the Hamiltonian. If we turn off all the 
perturbations, namely m, /i, h+ and then there are 
four such candidates falling into the second category. To- 
gether they close a C/(3) x U(l) algebra. The three mu- 
tually anti-commuting matrices, closing the CI (3) sub- 
algcbra act like standard two dimensional Pauli matrices. 
The remaining one, belonging to the U(l) commutes with 
three matrices, which close the C7(3) sub-algebra. With 



underlying topological pairing the CI (3) x U(l) algebra 
is constituted by 

M T = {to <g> 0:3, t ® j3,r ® r, r <g> iaia 2 } , (74) 

where the last entry belongs to the U(l) part. How- 
ever, due to Nambu's particle-hole doubling of the origi- 
nal problem, an 8-dimensional fc-dependent perturbation, 
afeMfe can acquire a finite expectation value, only if the 
matrix satisfies the algebraic constraint 

M k = T (n ® h) (n ® h) . (75) 

In the above equation, the =p signs depend on whether 
the coefficient cik is even or odd under the parity 
transformation (k — ^ — k). If the coefficient Ofc is of even 
parity, the only matrix satisfying the above constraint is 
To <£> r. On the other hand, for an odd parity (e.g., 
linear in k) the remaining three matrices can acquire 
finite expectation values from the zero energy sub-space. 
Notice that one of the matrices, t ® 0:3, appearing in 
the C7(3) part of Mt, multiplies k z in H gen , shown in 
Eq. (|60p . Therefore, the momentum along the vortex 
core k z does not cause any mixing of the zero energy 
states with the rest of the spectrum. Consequently, a 
first order perturbation calculation in terms of k z leads 
to the exact solution of the dispersive modes along the 
vortex core. The matrix To <8> CK3 acts as an off-diagonal 
Pauli matrix in the zero energy subspace, and hence the 
one dimensional dispersive modes are the symmetric and 
the anti-symmetric combinations of the fermionic zero 
modes with underlying point vortex. The energies of 
these two dispersive modes are E = ±fc z . When k z = 0, 
we have two Majorana fermions, which hybridize via k z 
and become complex fermions. 

The exactness of the perturbation theory in terms 
of k z is also applicable when we take into account 
the perturbations m,fJ,,h±. However, not all the ma- 
trices in Mt anti-commute with the generic Hamilto- 
nian H gen [0, At'jHz = 0]. Before, we proceed to prove 
this statement it is worth appreciating an algebraic 
identity. 29 ' 36 ' 64 Expectation value of an operator (M. ) can 
be expressed as 

^>4( e - e y^E, (76) 

^ occupied empty ' 

where \1/ e are the eigen-states of a generic Hamiltonian, 
%. If there exists a matrix, say T, which anti-commutcs 
with 7-L and commutes with Ai, the above mentioned 
sum is restricted to the zero energy subspace. When 
m = /x = h± = 0, and M. = (tq ® as)k z , one can 
choose T = To (8 iaia 2 - When m,fx,h± are finite, we 
can still find a matrix, T with requisite criteria, for 
the chosen M = (tq x a^)k z . If all the perturbations 
to, (i, h± are nonzero, we cannot find any unitary matrix 
for T. Rather there is an anti-unitary operator (At), 
namely (ti <S> ifiaxa?) K, where K is the complex 
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conjugation, which anti-commutes with the Hamiltonian 
H gen [0, At] k z — 0] and commutes with (t ® «3)fc 2 - 
Therefore, we can choose T = At- Hence, the one 
dimensional dispersive modes are always the symmetric 
and the anti-symmetric combinations of the zero energy 
modes bound to the point vortex. Even in the presence 
of the gauge holds, we can still choose T — At, and the 
above conclusions remain unaltered. 



B. s-wave pairing 



constructed from the Majorana zero modes due to the 
point vortex. Here, we first present the exact solutions 
of the dispersive modes for line-vortex as well as the 
zero mode solutions of the point vortex, for particular 
choices of the parameters in H gen [A s , At]. Then we 
show that the one dimensional dispersive modes with 
energies E = ±fc z , are either linear combinations of 
(for the topological pairing) or exactly (for the s-wave 
pairing) the zero modes for the point vortex. Some 
particular limits of this problem has been considered 
previously in Refs. Il6l - fl9l . 



A similar conclusion can be arrived at even with an un- 
derlying s-wave pairing if we turn off all the chiral symme- 
try breaking perturbations, for example to, Further- 
more, when /i, h + is set to zero the C7(3) <g> U(l) algebra 
of the matrices, anti-commuting with the Hamiltonian 
H gen [A s , 0; k z = 0] is constituted by 



{T Q (g>T, r (g) a 3 ,r 3 ® /3,t ® iaia 2 } , (77) 



where the last entry belongs to the U(l) part. Ap- 
pearance of the matrix To <S> 03 in the C7(3) part of 
M s , immediately guarantees that the dispersive one 
dimensional modes can be obtained by performing a 
perturbative calculation over k z within the two di- 
mensional basis spanned by the fermionic zero energy 
modes due to a point vortex in the k z = plane. It 
can also be confirmed from the Eq. (|76[) , upon choosing 
T = tq (8 ia\OL2- Otherwise, the matrix to ® «3 acts 
as the diagonal Pauli matrix. Hence, the dispersive 
modes, with energies E = ±k z , arc identical to the 
fermionic zero mode due to the point vortex. Let us 
now incorporate a finite chemical potential (fx) and 
the symmetric Zeeman coupling h+. One can then 
choose T = (T2 £8 i/3a\a2) K, for M = (t (g> az)k z . In 
conjunction with such choice of T, Eq. (|76|) guarantees 
that the dispersive modes with E = ±k z , are exactly the 
two fcrmoinic zero modes bound to the point vortex. 

It is worth mentioning that one of matrices in the 
C7(3) subgroup, namely T3 <S> /3 appears in H gen in Eq. 
(|60|) with the Dirac mass (to). This matrix also satisfies 
the condition in Eq. (|75[) . if its coefficient is momentum 
independent. Therefore, the internal structure of the 
zero energy sub-space shows that the Dirac mass is 
sufficient to cause splitting of these two states and place 
them at finite energies, ±to. 



V. EXACT AND PERTURBATIVE SOLUTIONS 
OF LINE VORTEX 

The internal structure of the zero energy modes of the 
point vortex allowed us to show that the one dimensional 
dispersive modes along the core of the vortex can be 



A. Topological pairing with /i = h+ = 

Let us first present the solution of the dispersive one 
dimensional modes along the vortex core with underlying 
topological pairing, when fi = h+ = 0. The coupled 
differential equations for gaplcss modes along the vortex 
core read as 

1 

— 1 

r 



(78) 



(to - /i_)Aj + (-i)e^ d r + -84, - a J A+ 



(to + /i_)A+ + (-i)e-** d r - -84, + a J A 



+ fc z Af+A r e-^ A; =EA+, 



k z A+-A r e- i *(A+) =EA-, (79) 



(m - h-)A^ + (-i)e # (d r + -8$ - a^j 



fc z A--A r e~^ A~ =EA+, (80) 



(to + fc_)Af + (-i)e-** (d r - l -d^ + a^j A+ 



fc,Af + A r e-** A+ =EAl 



(81) 



The remaining four equations can readily be found by 
taking the Hermitian conjugate of the above four. A r 
captures the radial variation of the superconducting gap, 
which in the presence of line-vortex vanishes smoothly 
as r —> 0, and saturates at A (constant) as r —> 00, 
otherwise arbitrary. Next we wish to find the solution 
for two different cases, i) E = k z and ii) E = —k z . In 
the former situation, the condition E = k z leads to a 
set of constraints among the various components of the 
Nambu-Dirac spinor, namely 



A: 



A+,A l 



-A+(Ar)t = (A+)t I (A7)t 



(82) 
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Upon imposing these constraints the above four equa- 
tions reduce to only two, and they read as 

{-i)e~^(d r ~ -dj, + a^AV + {m + h-)A+ 



A r e-^A^ =0, (83) 



{-i)e t +[d r + -d 4> 



At - (m - h-)A~ 
t 



A r e" 



A 



0. 



(84) 



These two coupled equations can be solved using the 
ansatz 

A+ = K(r) e-** g(r); A7 = TV (r) /(r), (85) 

where TZ(r) = exp (—if — Jq A r >dr'). We here take the 
profile of the vector potential a to be same as in Sec. II 
D. With such a profile of a, one obtains 



g(r) = a h 
f(r) 



t\J m 2 — h 2 

\J m 2 — h 2 



m + h 
Ci\l r Iq 



m — h 

when r<(. Outside the core of the vortex (r > £) 



(86) 



<?0) 


= c 3 I 1/2 


r\m\ 


+ cj. 


-1/2 


r\m\ 


m 


= ah/2 


r\m\ 


+ c 3 I_ 


-1/2 


r\m\ 



(87) 



The solutions and their first derivatives need to be 
matched at r — £, where the solutions for r < £ can be 
found by replacing h 2 — > h 2 + 1/2 A 2 4£ It eliminates two 
out of three arbitrary constants from f(r) and g(r), while 
the remaining one is fixed by the normalization condi- 
tion. The complete solution in the basis of 8-componcnt 
Nambu-Dirac spinor shown in Eq. (|58[) reads as 



\E 



if +ik z z - Jg A r ,dr' 



.(88) 



/ g(r)e-* \ 
if(r) 
-g(r)e-^ 
if(r) 

f(r) 
-ig{r)e l * 

V m J 



On the other hand, when we wish to solve the dispersive 
mode with E = —k z , the components of Nambu spinor 
are related according to 

Af = -A+ AT/ = A+ (Apt = _(A+)t,(A7/)t = (A+)t. 

(89) 



These constraints also reduce the number of coupled 
equations from four to two. Otherwise, the solution in 
terms of A^ and A7 are identical to that for E = +k z 
mode. In terms of 8-component Nambu-Dirac fermions 
the dispersive gapless mode with E = —k z assumes the 
form 



/ ff(r)e-* \ 
-if(r) 
g(r)e-^ 

if(r) 
ig(r)e t< >' 

-m 

ig(r)e lrj > 
V f(r) J 



\E 



i+ik,z - /„*■ A r ,dr' 



(90) 



where the radial functions f(r) and g{r) are identical to 
the ones in Eq. (JSBJ and Eq. ([87]). for r < £ and r > £, 
respectively. 



Next we present the solution of the fermionic zero 
modes with an underlying point vortex, when fj, = h+ = 
0. With this particular choice of parameters, one set of 
coupled differential equations for the zero mode reads as 



'+ (d r --d 4> + a<j> A7 + i'(m + /i_)A+ 



+ ^A r e- ^ *(A7) 



(91) 



ie^fdr + fa - <^Ja| + (m-/i_)A7 
+ A r e-*+(Af') =0. 



(92) 



The fermionic zero obtained by solving this set of equa- 
tions, in the basis of 8-component Nambu-Dirac spinor 
reads as 



/ g(r)e-** \ 



■ J r A r ,dr' if{r) 
ig(r)e^ 



V f(r) J 



e «e 



(93) 



where the two radial functions f(r) and g(r) takes the 
identical form as in Eq. and Eq. ([57]) for r <C £ and 
r > £, respectively. The other set of coupled differential 
equations for the zero mode reads as 



-itp 



d r - -d<j, + J Aj" - i(m + h-)A^ 

t 



iA r 



A; 



= o, 



(94) 



13 



A r e~ 



A: 



(95) 



The corresponding zero energy state takes the form 



/ \ 

-if(r) 




-/« 
ig(r)e^ 




(96) 



where f(r) and g(r) are identical as for the other zero 
mode, |*?). 

It is now evident that the solution for two one di- 
mensional dispersive modes, \E = +k z ) and \E = — k z ) 
in Eqs. (|88jl and (|90|) . are respectively the symmetric 
and the anti-symmetric combination of two fermion zero 
mode in Eqs. (|93f and (|96)l , in the presence of point 
vortex. 



B. s-wave paring with m = fi — h± — 

Let us now consider a line vortex along the z-direction 
with underlying s-wave paring. The dispersive mode can 
be solved analytically when we set m = fx = h± = 0. 
Furthermore, we also turn off the orbital contribution of 
the gauge field, i.e. a<j, = 0. Then the coupled differential 
equations read as 

i 
r 
i 

—i 

r 

(-*) (d r - ±d+) A+ + k z A+ + e-^Ar \AJ ) = £A f , 

(-*) (d r + l -d^j A+ - k z A+ + e-»Ar (Af ) f = EA7\ 

(97) 

The remaining four equations are the Hermitian conju- 
gate of above four. Upon imposing the constraints over 
the spinor components, as shown in Eq. (|82j) . the above 
four equation reduces to two decoupled equations. For 
E = +k z they read as 

ie~^ fd r - ~dA Al + e-^Ar (a;) T = 0, (98) 



(-*) ( d r - -d+ ) Al + fc 2 Af + e -*Ar (A+ j = EA+, 

+ )=EA+, 
t 



ie*+ ( d r + -<9 A+ + e"^A r (A 



0. (99) 



The first equation yields 



whereas the second one gives 



C 



A+= r± e -^e-So r ^dr' 



(100) 



(101) 



where C± are the arbitrary constants. However, to keep 
the second solution well behaved near the origin, we have 
to set = 0, therefore A^ = 0. On the other hand, the 
constraints in Eq. ([81?)) yield the following two decoupled 
equations for E = —k z 



(-i)e-^(d r - l -d^A- i+ e-^A r (Aiy = 0, 



(102) 



(_i) e <* fd r + -dA A+ + e-^Ar (a+) T = 0, 



giving 



and 



A7 = (i) CZe-^e-K^'dr'^ 



a 



A+ = —t e i %e-f° A r' dr ', 



(103) 



(104) 



(105) 



where arc also arbitrary constants. For the solutions 
to be well behaved in the vicinity of the origin we need 
to set C'l = 0. In conjunction with the constraints, the 
above solutions in Eq. (|100p and Eq. f|104[) . yield two 
dispersive modes with underlying s-wave order 



\ + k z ) = C±1Z{r,z) 



,\-k z ) = ClK(r, z) 



( \ 

-1 


1 


—i 


V i J 



where 7Z(r, z) = exp (— i-j — A r idr' + ik z z). 



i 


i. 



1 



v i y 



(106) 



Next we proceed to obtain the zero energy (E = 0) 
modes with an underlying point vortex of s-wave par- 
ing in the k z = plane. It can readily be solved from 
Eq. (J97J upon setting E = 0, k z = 0. In light of the 
above discussion, readers can convince themselves that 
Aijt", A^ and their conjugates in the zero energy subspace 
are proportional to 1/r. Therefore, these components are 
exactly zero for the normalizable zero mode solutions. 
The remaining components are finite for two zero energy 
modes. For one Majorana mode 



AT 



A 



(107) 
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and, for the other one 

A7 = -A+ = CZe- i *e-& A *" ir '. (108) 

k z acts like 173 = diag.(l, —1) matrix in the zero energy 
sub-space. Therefore, it does not cause any mixing be- 
tween two Majorana modes. Consequently, the gaplcss 
modes are same as the vortex zero modes, multiplied by 
the plane wave factor exp ik z z. 



VI. DISCUSSIONS 

In this paper we have demonstrated that the original 
Jackiw-Rossi Hamiltonian describing a point vortex 
defect in two spatial dimensions (see Eq. (JTJ), can be 
augmented by two additional terms (see Eq. (JSJ), which 
still possesses a spectral symmetry. Consequently, the 
generalized Jackiw-Rossi Hamiltonian of Eq. ([5J may 
support zero modes. For generic values of the pertur- 
bation parameters satisfying the condition in Eq. (|48|) . 
we obtain a single zero mode only for the odd vorticity. 
In contrast, there are no zero modes for the even 
vorticity, for generic perturbations. To demonstrate the 
emergence of the Z2 index for the zero modes, we have 
employed a method of matching asymptotic solutions, 
which correctly captures all the known results within a 
single framework. We have also found that there exist 
special values A = ±x, for which there are n number of 
zero modes respectively for a vortex and an antivortex 
of vorticity n. 

One of the main goal of this paper is to determine the 
number of gapless one dimensional modes along the line 
vortex of a gapped paired states of three dimensional 
Dirac quasiparticles. In order to answer this question, 
we have mapped the problem into the determination of 
the number of vortex zero modes of appropriate iJ^j's. 
Through this procedure, we have succeeded in showing 
that the number of gaplcss modes is also dictated by 
the Z 2 index of if^'s. We have exemplified this by 
considering the topological odd parity and trivial s-wave 
pairings. If the underlying Dirac fermions are massless, 
then in the absence of the chemical potential and the 
Zeeman couplings, both types of pairing lead to two 
copies of appropriate Hjr for k z = 0. Consequently, 
the number of the gapless modes is governed by the Z 
index theorem of Weinberg. When generic perturbations 
are considered, only the topological pairing sustains two 
gaplcss modes for odd vorticity, in accordance with the 
Z2 index. 

Our calculations for isolated vortex can be extended 
in perturbation theory set up for multiple vortices in 
the dilute vortex limit. In the presence of multiple 
vortices, if we again consider k z = 0, there will be 
tunneling within each copy of H^^'s (see Eq. (1551) and 



Eq. (|67|) 'l for topological superconductor. As far as 
the topology of the order parameter field is concerned, 
there is no difference between a n-vortex and widely 
separated n-number of single vortices. Our analysis now 
suggests an interesting even- odd effect based on the 
Z2 index. When the net vorticity is odd, the gaplcss 
state will survive the tunneling effects. In contrast, the 
gaplcss state will be absent for a net even vorticity. This 
consideration can also be extended to non-dilute limit 
following the strategy in Ref. |65|. 

Our analysis can be applied to an interesting problem 
of chiral anomaly for the gapless one dimensional modes 
along the vortex line of an axial superfluid, which was 
considered by Callan & Harvey^ They have considered 
the following model 

3 

Hax = ^2 7o7j (-idj - eAj) + A(r)7 exp (16*7075) , 

. (109) 

where we have used five mutually anti-commuting 7- 
matrices, and Aj is the electromagnetic vector poten- 
tial, and e is the electron's charge. When, we consider 
a line n-vortcx of the axial superfluid order parameter 
along the z direction {9 = n<f>), there are n gapless chiral 
one dimensional modes. This number of modes is tied 
with the Weinberg's Z index theorem for the underly- 
ing Jackiw-Rossi problem in the x — y plane. The chiral 
one dimensional modes in the presence of the electric 
field along the z direction, give rise to a non-dissipative 
electric current along the z direction, determined by one 
dimensional chiral anomaly j z — n x e 2 E z /(2ir). This 
current in turn is supplied radially from the bulk into 
the vortex core, which is captured by the following axion 
electrodynamics term 

C axml =--^ J d A xe^P x d^9 A v d p A x . (110) 

Now we may add various fermion bilinears to the above 
model, which can still support gapless modes along the 
vortex line. According to the construction of in 
the previous sections, we can add 75 A and 17172X1 which 
respectively describe an axial chemical potential, and 
the third component of the space-like axial vector (these 
terms break the Lorentz and CPT symmetries). Now 
there is a single gaplcss mode only for the odd vorticity, 
if Ag + A 2 > x 2 ■ The value of the amplitude A at radial 
infinity has been chosen to be A . Consequently, we can 
find non-dissipative current only for the odd vorticity, 
under generic values of these parameters. Accordingly 
the bulk axion term, which is usually computed through 
the Goldstone-Wilczek formula,— has to be modified to 
capture the Z2 index. 

The condition A 2 , + A 2 > x 2 has a simple physical 
meaning in terms of the uniform model's band structure. 
If x — 0; the Kramer's degeneracy is lifted by A, but 
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keeping the spectrum fully gapped. On the other hand, 
for A = 0, the spectrum is fully gapped only when 
A > x- F° r X > A we have a Weyl semi-metal phase, 
which does not support the vortex zero modes. In 
the Weyl semi-metal phase, there are chiral surface 
states, which lead to anomalous Hall conductivity and 
chiral magnetic conductivity^ The anomalous transport 
properties in the gapless phase are also captured by 
appropriate axion electrodynamics terms, which are not 
related to vortex zero modes. This is also interesting 
to note that the number of gapless modes also controls 
an associated gravitational anomaly formula, and our 
work suggests its Z 2 modification in the presence of A 
and/or \- We also note that the axial vector % breaks 
the spatial rotational symmetry, and for this reason the 
zero modes can only be found for a line vortex aligned 
with the axial vector. 
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Appendix A: Zero energy states at A = ±x 

We here present the detail solutions of the zero energy 
states with an underlying anti- vortex for A = ±%. When 
A = — x, one can write the the equation for the zero 
energy modes for even-vorticity (n = 2s) from Eq. (fTB])- 
(TlTJl) . Written slightly differently, they read as 



d r + ^-jf 1 (r)+A r g 1 (r)=0, 
d r - ^jf 2 (r)+A r g 2 (r) - 2A.fr (r) = 0, 

- 2s ~ 1 ~ p )g 2 ( r ) + A r f 2 (r) - 2X 9l (r) = 0, 



zero energy modes can only be found when c\ — c 2 = 0. 
On the other hand, in the vicinity of the origin 



d r + ^— ^ ) gi (r) + A r /i(r) = 0, 



(Al) 

where p is a positive definite integers, and takes the values 

p = 0,l,..-,2«-l. (A2) 

In the vicinity of the origin, where A r —> 0, ff~(r) = 
&ir~( p+1 \ and gf{r) — c 2 r~( 2s ~ p \ Hence, normalizable 



(A3) 



/2 W = c i rV , g 2 ( r ) = c ir 

are well behaved functions for all p given in Eq. 
Far away from the origin, these two functions are f 2 (r) = 
9 2 ( r ) = c ex P (~ A ?') . Two out three arbitrary constants 
(ci, C2, c), can be fixed by imposing the boundary condi- 
tions 

fi(r = = /f> = 0. 9i(r = = 9?(r = 0, (A4) 

while the remaining one is determined by the overall 
normalization factor. Hence, there are n = 2s number 
of zero energy states with underlying anti-vortex, when 
A = — x- In terms of the original functions in Eq. Q14[), 
V(x) = for the zero energy modes. Therefore, the zero 
energy modes are eigenstate of the operator Ts = <jq ® 0-3 
with eigenvalue +1. These solutions match exactly with 
the ones found by Jackiw and Rossi, in the absence of 
gauge fields 3 

If, on the other hand, A = x, we have 



/<(r) = d 2 r p , /<(r) 



c 2 A 
p+l' 



g<(r) = c' 2 r 2 ^-^ g<{r) = 



p- 2s 



r 2s - p , (A5) 



where c 2 ,c 2 are arbitrary constants. Far aways from the 
origin the functions behave as 

f 2 { r )=9 2 { r ) = cexp(-A r), 

fi(r) = g?(r) = -2cXr exp (-A r), (A6) 

where c is also an arbitrary constant. After satisfying 
the boundary conditions, for example the one shown in 
Eq. (|A4[) . two out of three arbitrary constants, (c 2 , c' 2 , c) 
are fixed. With only one arbitrary constant, it is now 
impossible to satisfy two similar boundary conditions 
for f 2 (r) and g 2 (r). Therefore, when A = Xi there is no 
zero energy state when the underlying anti-vortex has 
even vorticity. 

Let us now consider the anti- vortex with odd vorticity, 
n = 2s + 1. We first focus on 2s ansatz of the form Eq. 
(|14[) . The equation of the zero energy modes then reads 
as 



+ ^|/ 1 (r)+A,g 1 (r)=0, 



/ 2 (r) + A r52 (r)-2A/i(r)=0, 



d r - ^r^).92(r) + A r / 2 (r) - 2A ffl (r) = 0, 
+■ 2S ~ P+1 \ gi (r) + A r .h(r) = 0, (A7) 
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when A = — %. For the normalizablc zero energy modes, 
we find /i(r) = gi(r) = 0. The remaining two functions 
in the vicinity of the origin read as 

f<(r) = Cl rP, g<(r) = c 2 r 2s -P, (A8) 

and far away from origin they are 

fit?) = 5 2 > W = cexp(-A r). (A9) 

Hence, for A = — x, there exists 2s number of zero energy 
states with an underlying anti-vortex of vorticity 2s + 1 
of the form Eq. (| 14|) . These solutions are also identical 
to the ones one find when A = x — 0. 20, However, when 
A = %, there is no zero energy mode of the form Eq. (|14[) . 

Next we solve the zero energy mode with the ansatz of 
the form Eq. (|39|) . when A = ±x- One can also write the 
Eq. (|4Tj) for the zero energy states as 

+ + A -/M + ( A + *M r ) = o, 

[d r - ^ f(r) + A r f(r) - (A - x )g(r) = 0, (A10) 



where p = (n — l)/2, and n(odd) is the vorticity. For 
A = —x, the radial functions of the zero modes are given 
by 

f(r) = 0, g(r) = cr"exp^-| A,„ dr'^j , (All) 

where c is an arbitrary constant, similar to the one found 
in original work by Jackiw-Rossi.— With this particular 
ansatz, we also find normalizablc zero energy modes even 
when A = x- The radial functions then read as 

f(r) = -c (^y) r p+1 cxp (-£ A r > dr') , 
g(r) = cexp|^-| A r , dr'^j . (A12) 

Therefore, the zero energy mode of the form Eq. (|39p . 
exists whether A = x or ^ — ~X- F° r n = 1 or p = 0, 
this solution matches exactly with the one shown in Ref. 

a. 
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